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Recent experiments have demonstrated an open system realization of the Dicke superradiance 
transition in the motional degrees of freedom of an optically driven Bose-Einstein condensate in a 
cavity. Relevant collective excitations of this light-matter system are polaritonic in nature, allowing 
access to the quantum critical behavior of the Dicke model through light leaking out of the cavity. 
This opens the path to using photodetection based quantum optical techniques to study the dy- 
namics and excitations of this elementary quantum critical system. We first discuss the photon flux 
observed at the cavity face and flnd that it displays a different scaling law near criticality than that 
obtained from the mean fleld theory for the equivalent closed system. Next, we study the second 
order correlation measurements of photons leaking out of the cavity. Finally, we discuss a modu- 
lation technique that directly captures the softening of polaritonic excitations. Our analysis takes 
into account the effect of the finite size of the system which may result in an effective symmetry 
breaking term. 
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I. INTRODUCTION 

In the study of strongly correlated systems and collec- 
tive phenomena, light has traditionally assumed the role 
of a spectroscopic probe. Recent progress in the con- 
trol of light-matter interactions through cavity QED has 
brought forth new systems where light and matter play 
equally important roles in emergent phenomena. Such 
hybrid light-matter systems are characterized by the ex- 
istence of well-defined quasi-particles, polaritons, which 
are partly light partly matter-like. A crucial feature of 
these systems is that they are inherently out of equi- 
librium due to unavoidable photon-leakage, giving rise 
to open system analogues of certain well-studied quan- 
tum many body Hamiltonians. Such systems may for 
instance be formed by scaling up standard single-cavity 
QED systems to lattices of cavities. Recent theoretical 
work in Cavity QED lattices has addressed the realization 
of superfluid-Mott insulator transition of polaritons [IHl] , 
fractional quantum Hall states [S], the Tonks-Girardeau 
gas in one-dimensional geometries [51 [7] and effective spin 
models [3l [8] . 

A different venue is to couple non-trivial states of mat- 
ter to cavities, where the formation of polaritonic col- 
lective excitations may result in new emergent phenom- 
ena. One particular approach is to couple a Bose-Einstein 
Condensate (BEC) to a single mode of a high-finesse op- 
tical cavity. This results in tunable, long-range forces 
between atoms of the BEC that is mediated by the cav- 
ity field. A phase transition from a uniform BEC to 
a self-organized, density-modulated phase has been pre- 
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dieted 151 as a function of the power of the laser driving 
the atoms transverse to the cavity axis and experimen- 
tally observed recently [TD]. A similar self-organization 
transition has also been observed for a thermal cloud of 
atoms [llj . It remains a question to distinguish the two 
types of transitions from each other, a matter that can 
be settled by studying the nature of excitations around 
criticality. 

Here we study the excitations of the zero-temperature 
BEC self-organization transition, which, in the experi- 
mentally realized parameter regime, was shown to be a 
faithful open-system realization of the Dicke superradi- 
ance transition where spins correspond to the collective 
motional degrees of freedom of atoms. The transition 
is driven by softening of a polaritonic excitation as the 
critical point is approached. This provides access to the 
internal excitation dynamics close to criticality through 
photons that leak out of the cavity mirrors. Thus well- 
established quantum optical measurement schemes ap- 
pear as specially suited to monitor the intra-cavity exci- 
tation dynamics as well as the phase diagram. We first 
discuss the critical behavior of photon flux measured at 
the cavity mirror. We then study the signature of excita- 
tions in the second order correlation functions of leaking 
photons, a measurement that can readily be performed 
using standard quantum optical schemes such as the 
Hanbury-Brown-Twiss setup [12] . Finally, we consider 
a modulation technique that directly captures the soft- 
ening of the relevant polaritonic mode through the pho- 
todetection of leaking photons. Our discussion will take 
into account the role of the finite size of the system that 
is relevant for a realistic experimental setting which is 
shown to result in an effective symmetry-breaking term. 

The paper is organized as follows. In Section |ll] 
we summarize the system and the governing Hamil- 



tonian together with mapping to the Dicke model. 
We then investigate the collective light-matter excita- 
tions in Section |III[ Finally we discuss three differ- 
ent photodetection-based measurement schemes in Sec- 
tion HVl 



II. MODEL OF THE SYSTEM 

We consider a Bose- Einstein condensate (BEC) of 
length D coupled to a single mode of frequency Wc and de- 
cay rate k of a high finesse cavity of length L, driven by a 
laser with frequency tOp from a direction perpendicular to 
the cavity axis (TUl. We assume that |Ac| = |ajp — wd ^ k. 
so that the cavity can be quasi-resonantly excited by the 
scattered pump photons from atoms. Simultaneously, the 
laser is red-detuned far from an internal atomic transi- 
tion at uJa, so that |Aa| = \ujp — U!a\ 3> 7a, where 7^ 
is the atomic linewidth. This ensures that the atoms 
are predominantly in their ground states during the ex- 
citation process, suppressing spontaneous emission and 
giving rise to an optical potential for the motional de- 
grees of freedom of atoms \E'^{x)\'^/Aa. Here E~^{x) — 
gQ(Pc{x)a+flp(pp{x) is the positive rotating component of 
the electric field felt by an atom of the cloud at position 
X, due to the interference of the cavity field (photon cre- 
ation operator a, atom-field coupling gg, mode function 
(pc{x)) and a coherent laser field (with Rabi frequency 
rip and standing wave pattern ipp{x)). The Hamiltonian 
under these approximations is given by 



H = -hAca'^a+ / dxi'\x) 
+ V{x)]i'{x) 



\E+{, 



2m 



Aa 



(1) 



We consider a situation where the BEC is trapped by 
an additional external trapping potential V{x) and that 
this is deep in the radial direction confining the cloud 
along the cavity axis (x). With the additional assump- 
tion that the driving laser beam is broad, we can reduce 
the problem to a an effective one-dimensional problem 
with (pp{x) w const, (the constant to be absorbed into 
rip) and (pc{x) w ipc{x) [l3]. We will assume that the 
cavity mode function is given by (pc{x) = -jj sin(Ga;) 



where G — I-k jX^ 
then given by 



^■n jXp. The final Hamiltonian is 



H = -hAca}a+ \ dx^^(x) 



->U 



-fi2 d2 



V{x) 



2m dx^ 

+ hUo\(pc{x)\'^a'^a + hrnpc{x){a^ + a)] *(x 



(2) 



where C/g = 9o/Aa and 77 = flpgo/Aa- We have sub- 
tracted the energy provided by the constant potential 
hVti/Aa- We will neglect the contact interactions be- 
tween atoms of the condensate, which is not essential for 
the physics discussed here. 

As a function of the tunable pump power 77, this model 
displays a phase transition from a non-organized phase. 



a homogenous condensate, to an organized, density- 
modulated phase O [m [15] . Below a threshold power 
r^c, the intracavity (mean) field is vanishingly small and 
the cloud is in the ground state of only the external trap- 
ping potential. As the critical point is crossed, the atoms 
self-organize into a crystalline order. This in turn results 
in a non-zero cavity mean-field through the scattering of 
the pump photons from the density-modulated atomic 
cloud into the cavity. In the organized phase the system 
chooses spontaneously between two density-modulated 
configurations that are shifted by half a cavity wave- 
length: the even or odd checkerboard configuration. This 
self-organization phase transition was observed experi- 
mentally for thermal atomic gas [11] and for a cloud of 
atomic BEC [T^. 

For the experimental conditions of Ref. [j^, the self- 
organization transition can also be seen as the Dicke tran- 
sition from a normal to a superradiant phase [10| , ll6 j . We 
will confine ourselves to this regime but take steps to ac- 
curately model the experimental conditions. We assume 
a BEC considerably smaller than the cavity size, imposed 
by the external trapping potential V(x) and expand the 
field operator 



^{x) 



i4>n{x) 



(3) 



Here, (/>„(a;) is the atomic single-particle basis satisfy- 



ing 



2m da;' 



+ V{x] 



We will as- 



4>n{x) = UJn4>n{x). 

sume V{x) to impose Neumann Boundary conditions at 
X = {L - D) /2 ± d a.nA X ^ {L + D) /2 ± d, allowing for 
an asymmetric placement of the trap by a length d with 
respect to the cavity walls (we will assume d/D ^1). 
Then, (j)o{x) = l/y/D is the uniform mode of the conden- 
sate with zero momentum. In the expansion p|, we keep 
only one additional mode having a relatively large over- 
lap with cavity mode ipc{x), say 4>n{x), which becomes 
the dominant contribution in the Hamiltonian of Eq. ([2]) 
when the cavity mean-field is non-zero. The correspond- 
ing wave-vector is fc„ ~ G. Here it is assumed that these 
two modes satisfy the relation CgCo-l-cJjC„ = N. Introduc- 
ing the Schwinger representation for 'spins' J_ = CQCn, 



J+ = J_ and Jz = (cj,c„ — cJco)/2, the final Hamiltonian 
can be written as 



uid* a 



+ ojoJz + -^{a'^ + a){J+ + J^) 



A 



+ -^(d^ + d) 



N 



J, 



(4) 



where uj = —Ac + {NUq/ D) J dx\(fc{x)\'^ , ujq = ujb. = 
hG^/2m, A = iy/ N/D)r] J dxipc{x)(j)nix) and A' = 
{y/N/D)r] J dxipcix)4>oix)- In all of our calculations 
here, the limits of the overlap integrals are those of the 
atomic trap since (j)n{x) = outside. Since we are in- 
terested in the critical behavior, an additional dispersive 
shift term |10j which is negligible near critical point is 
dropped in Eq. Q. 



For A' = 0, Eq. [4] is the well-known single-mode Dicke 
Hamiltonian. In the thermodynamic limit of iV 3> 1, the 
Dicke Hamiltonian m exhibits a quantum phase transi- 
tion at a critical coupling strength Ac = (l/2)^a;wo from 
a normal phase with (a) — (J-) = to a superradi- 
ant phase fi7\^^ with (a) 7^ 0, (J_) ^ 0. The parity 
symmetry of the Dicke Hamiltonian (|4| under (a —>■ —a, 
J_ — > —J-) gives rise to two distinct degenerate super- 
radiant phases that break the parity symmetry of the 
Hamiltonian. 

We observe that a non-zero A' acts as a symmetry- 
breaking bias field where the sign of A' determines which 
state is chosen in the superradiant phase. If a trap is not 
positioned perfectly symmetrically with respect to the 
middle of the cavity, A' ^ 0. This was observed in recent 
experiments studying the process of symmetry breaking 
in real time through an interferometric heterodyne de- 
tection scheme [25] . 

Finally, cavity losses at a rate k can be included 
through a Lindblad master equation approach |26j . 
Note that this leads to an open-system analogue of 
the Dicke superradiance transition, which, most signif- 
icantly results in a shift of the critical point to Ac = 
(l/2)-\/(a;o/w)(K^ +aj^). A second important conse- 
quence is that the excitation spectrum and the dynamics 
becomes dissipative. 

Introducing mean fields a — (a), /3 = (J_) and w — 
{Jz) in the Heisenberg equations of motion for the system 
Hamiltonian Q, one gets 



A A' 

a — —(K + iuj)a — i^^(B + /3*) — i^^ 



N 
~2 



(5a) 



threshold) disappears and the stable solution becomes 
nonzero for all A > values (Fig. [ljc,d)). We do not 
plot /3, it displays a similar behavior to a, with its sign 
opposite to that of a. Here, we should note that the 
spatial structure of the self-organized state, as shown in 
Fig. [IJb), depends on the sign of A'. 
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FIG. 1: (a) The real part of a for A' = 0. Below threshold the 
stable solution is a = and above threshold there are two sta- 
ble solutions with a phase difference n. (b) The steady state 
atomic density profiles for A' < and A' > with A = 9 and 
A'/A ~ ±120. (c) Dependence of a on the coupling parame- 
ter for the trap displaced to the right, which corresponds to 
sgn{X) — sgn{X'). (d) Same as (c) except that the displace- 
ment of the trap is to the left so that sgn{X) = —sgn{X'). In 
(a), (c) and (d) blue (red) lines represent stable (unstable) 
solutions. For all of the calculations, we use the parameters 
uj — 300, K = 200, ujo = 1 and for these choice of parameters 
Ac ^ 10.41. 



A A' 

B = -iuJoB + 2i^=(a + a*)w + i^=B(a + a*), (5b) 



III. POLARITONIC EXCITATIONS 



(a + a*)(/3-/3*). 



(5c) 



These equations have to be solved with the constraint 
that pseudo angular momentum |/3p -I- w^ = N"^ /A is 
conserved. Analytical solutions can be found for this 
set of nonlinear equations in the steady-state when A' — 
[26]. In that case, the steady state solution displays 
a bifurcation point at A = Ac. While ass = Pss = is 
the trivial solution for all values of A, it's only stable for 
A < Ac. For A > Ac, this solution becomes unstable and 
two new set of stable solutions appear given by [26| 



±^/N- 
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A4' 
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(6b) 



The behavior of these solutions for Ugs are shown in 
Fig. fT|a). For A' ^ 0, the bifurcation point (and the 



In this section, we investigate the collective light- 
matter excitations of the system that play an crucial role 
in the course of the underlying phase transition. To this 
end, we employ Holstein-Primakoff transformation to ex- 
press the atomic spin operators in terms of bosonic mode 

operators b and W such as J+ = wy N — Wh, J_ = jI_ 

and Jz = y^b- N/2 [Ml US]- Substituting these expres- 
sions into the Dicke Hamiltonian Q and expanding in 
the limit A^ ^ 1, we get the Hamiltonian governing the 
fluctuations around the steady state semiclassical solu- 
tions diss and Pss 

H„p/h = ujc^c + UJQ(Pd + gi {d^ + df + g2{c^ + c) {S + d) 

(7) 
where c and d are the photonic and atomic fluctuation 
operators respectively, a = ass + c and b = /3ss/VN + d, 
and 



lo'q = ujQ — 2A 



Ps 



1-^? 



=Re(ass), 



(8a) 



gi = —A = —iie(ass), 



.92 = A 



1-2/32 



1-/9.2 



-A'/3, 



(8b) 



(8c) 



Here we introduced the scaled variables ass = otss/vN, 
Pss = Pss/N . One should note that the steady-state val- 
ues ass and Pss are A'-dependent. The quadratic Hamil- 
tonian above leads to linear equations of motion h = Mh 
for fluctuations h = ((c), (c^), (d), {d))). This was done in 
the resonant case for oj — ujq and A' = in Ref . [26 . Here 
we are interested in the dispersive cavity regime of the 
self-organization problem T0','T6] for which Wq <C oj^-I-k^^ 
In that regime the real and imaginary parts of the lowest 
excitation eigenvalues are shown in Fig. [2] with A' = 0. 
The real parts of the remaining two eigenvalues are very 
large compared to those shown in the figure due to the 
dispersive nature of the cavity and are not shown. By 
employing perturbation theory for the small parameter 
e — lS^/{ijj^ + 't^), the real and imaginary parts of the 
polaritonic eigenvalue up to 0{£^) are calculated as 



woWl - 



A2 

a! 



1 



A^ 



2 w2 + k2 XI 



A2 



k2A2' 



(9) 
We observe that the energy gap monotonously decreases 
and the lowest energy mode "softens" as we approach the 
critical point from below. Note that this is the atomic 
excitation at A = which gradually acquires a photonic 
component as we approach the critical point, therefore 
we refer to this collective excitation as polaritonic. This 
is also the reason why the imaginary part of this excita- 
tion becomes larger, as in Fig. [2| as the critical point is 
approached. We note however that there is a very narrow 
regime around Ac in Fig. [2] where the excitation energy 
Heiujf^x) is zero and the damping lva{ujex) decreases to- 
wards the critical point. In this regime of critical slowing 
down fluctuations are overdamped. This regime is char- 
acterized by two values Ai < Ac and A2 > Ac below and 
above threshold respectively. For the dispersive cavity 
case, these two values can be approximated by 



Ai ~ Ac 



(w2 + ^2) 



2\2 



A9 — Ac 



l + :r- 



1 K^UJ, 



2, ,2 



2 (w2 + k2)2 

(10) 
The behavior of uj^x in this range is drastically different 
from Eq. (Bl) and is given to order (1 — \^ /X^) by 



mode provides a fluctuation channel that drives the non- 
equilibrium phase transition. Secondly, due to the in- 
creasing light-like content of this channel, we get a first- 
hand look into the fluctuations around the critical point 
by monitoring photons that leak out of the cavity. 
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FIG. 2: The real and imaginary parts of the eigenfrequencies 
of polaritonic excitations. Parameters are the same ones used 
in Fig. [1] The figures on the right column show the indicated 
magnified views of real and imaginary parts around A = Ac. 



IV. PHOTODETECTION MEASUREMENTS OF 
CAVITY PHOTONS 

Photon flux that leaks out of the cavity below the 
threshold is very low. Therefore, photodetcction mea- 
surements appear to be most useful for characterizing 
critical fluctuations. We consider below two measure- 
ment schemes, direct photon counting and second or- 
der photon correlations that can be implemented with 
standardly available photon counters. In section Sec- 
tion IV C we discuss a modulation scheme which pro- 



vides a direct access to the softening behavior of the po- 
laritonic mode. 

To this end, we make use of the standard input-output 
formalism j26| I27j , introducing the photonic input and 
output operators a^n, aout^ that couple to the intra-cavity 
fluctuation operators c and d in the following manner: 



(12a) 



-(iw -I- n)c — ig2{d + dJ) + v2Kdin, 



2n V A2 



(11) 



below threshold. 

The existence of this collective soft mode is signifi- 
cant for two reasons. Firstly, near the critical point this 



-iajgd - 



2ig2(d + S)-ig2{c + c^), (12b) 



with din being operator for the quantum noise incident 
on the semi-transparent cavity wall satisfying the com- 
mutations relation [ai„(Oi'^L(^')] ^ S{t—t'). For vacuum 
input, the relation (a^ (*)"««(*')) = {din{t)din{t')) = 



holds. The cavity output can then be expressed as 

aout{t) = V2K[c{t) + ass] - ain{t). 

A. Photon counting 

The photon flux measured outside the cavity can be 
expressed as Iss = {aiut{t)^out{t)) = 2k{c' {t)c{t)) . The 
photon flux was numerically found to diverge at the crit- 
ical point [TfiJ f26] but the precise form of scaling law was 
not elaborated on. We find that the intra-cavity photon 
number is given by 



{c\t)cit))ss = 



A2 



8A2 



1 - [X/XcY 



(13) 



clearly displaying non-equilibrium mean field critical 



scaling with an exponent j„ 



1. We note that this 



scaling is drastically different than that of the equilib- 
rium, ground state expectation value of the intra-cavity 
photon number: 



'/it 



cHt)m) 



A2 



fls 



a;2y^l-(A/A,)' 



(14) 



obtained through perturbation theory in e = uj'^/[uP + 
K^). The latter displays an exponent 7eg — 1/2 as ex- 
pected from standard equilibrium mean-field systems. 
The reason behind this is the depletion of the ground 
state (for T = 0) through coupling to the photonic en- 
vironment. This possibility was hinted on in Ref |16j . 
but the authors only calculated the rate of depletion in 
the short-time limit i.e. in the transient regime. We 
find here that the depletion settles at a steady state, giv- 
ing rise to an entirely different scaling law for incoher- 
ent fluctuations of intracavity photons as critical point 
is approached. Note that the two expressions do not 
even agree for k — >■ that enters the non-equilibrium ex- 
pression through the form of Ac. However, it should be 
pointed out that the limit k — >■ is a singular limit and 
has to be considered with care. The steady-state regime 
where this expression holds is shifted to i — ^ oo as k — >■ 0. 
This time scale can be calculated from the imaginary part 
of the polaritonic excitation branch and is found to scale 
as |Im(wea;)|^^ ~ uj'^ / {AujoX"^ k) according to Eq. ^. For 
K = 0, a steady state will never be reached unless the sys- 
tem starts out in the ground state. What remains to be 
seen is whether open quantum systems may be forming 
an entirely new universality class. 
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FIG. 3: (a) False color plot of the dependence of the steady- 
state second order correlation function Qsa (j) on A. We 
observe that the period of the oscillations increase with in- 
creasing A. Note that Qss \t = 0, A) — 3. (b) The loga- 
rithm of the Fourier transform of gL (color code) is plot- 
ted as a function of frequency i^ and the coupling parameter 
A. The two peaks other than the one at i^ = appear at 
1/ = ±2loo-sJI - (A2/Ai). In both figures A' = 0. The system 
parameters are same as in Fig. w\ (c) The long-time behavior 
of gls (t) displaying an exponential decay of the envelope of 
its oscillations. 



B. Second order correlations of photons 



The second order correlation function of leaking pho- 
tons g^"^^ (t, r) is given by 

(2)/^ ^ ^ {aluMalutit + T)dout{t + T)aout{t)) ,^^. 
(aLt(*)ao«t(t))2 



By using the input-output relations and the vacuum 
nature of the input noise, in the steady-state the expres- 



sion in Eq. ([T5|) can be simplified to 



5(2) (i,T) = 1 + l.g 



r)c{t)) 



„2 1 2 



2 1 a. 



((ct(t)c(i)) 



2\2 



,P) 



(16) 



where g'^'^\t,T) = {c\t + T)c{t)) / {c^ {t)c{t)) is the first 
order correlation function. This steady-state form can 



be calculated by solving the Eqs. (12 1 together with 
their adjoints in Fourier space. Some of the results for 
A' = are shown in Fig. Js] for different A < Ac val- 
ues. It can be seen that ^L (r) displays underdamped 
oscillations with the oscillation period progressively in- 
creasing as critical point is approached. A Fourier anal- 
ysis of the oscillations in Qss {t) as a function of A is 
shown in Fig. Islb). We observe that the position of the 
peaks Vpeak follow very closely the excitation frequency 
of the softening polaritonic mode given by Eq. (|9]); in- 
deed it can be shown that i'peafc(A) ~ 2Re[we2;(A)~ The 
width of the peaks, not resolvable for parameters chosen 
here, is proportional to Im[a;ea;(A)]. One can furthermore 
show that gss (0) = 3 as a consequence of the identity 
|(c'f(t)c(i))ss| = l(c(i)c(t))ss|, indicating that the evolu- 
tion of the system is into a two-mode squeezed state in 
the steady state. The above analysis provides a unique 
window into the nature of fluctuations close to the critical 
point. 

We next calculate gls {t) for nonzero values of A'. 
Here, the second order correlation function displays a 
beating pattern in time as seen in Fig. |4] This is the con- 
sequence of a non- vanishing mean-field Uss extending all 
the way below the threshold, providing a non-zero dc 
component ass in Eq. (16 1. Thus such a distinct beating 



pattern is a signature of the interference of a non-zero 
coherent cavity field and incoherent photons. 



C. Modulation Spectroscopy 

In this subsection, we analyze a modulation technique 
that provides an alternative access window into critical 
fluctuations. This technique relies on parametric reso- 
nances of the cavity-BEC system and is similar to mod- 
ulation techniques applied to the analysis of ultracold 
gases [28H3T] . 

While various parameters of the system may be mod- 
ulated, we choose one that appears to be experimen- 
tally most straightforward: modulation of the trans- 
verse pump power. We assume a periodic modulation 
of the pump Rabi frequency, ri{t) = rj[\ + ecQs{vt)] 
where v is the modulation frequency and e <C 1, re- 
sulting in the modulation of the coupling parameter 
A(i) = A[l + ecos(vt)] in the Dicke model, Eq. ^. Here 
we analyze the case A' = 0. The analysis below can 
straightforwardly be extended to non-zero A'. 

We assume k 3> Wq, so that the cavity field dynamics 
adiabatically follows the atomic dynamics and the photon 
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FIG. 4: Qsa {t) for A' — A/360 is represented by blue solid 
lines and for the case A' = (i.e. Uss = 0) by dashed red 
lines. From top to bottom, A = 2, 6, 8, 9, 10. The asymmetry 
in adjacent peaks is due to the presence of nonzero mean-field 
ass ■ The system parameters are as in Fig. [2] 



field a can be eliminated adiabatically. This results in the 
following equation for the atomic variable /3: 



f3 = -iujol3 + AiX'^{t)- 



1/2 



(^ + /3*), 



(17) 

where we used the stable solution for w which is negative. 
Writing (3 = f3ss + 5j3{t), assuming small fluctuations 
around the steady state, we obtain the following equation 

ioTu{t) = 5l3{t) + 5l3*{t) 

^u+[A-2e~cos(— t)]u = 0, (18) 

This is the well known Mathieu equation with A = 
1 — (A/Ac)^, e — (A/Ac)^e and i = uj^t. According to 
the Floquet theorem, the solutions to this equation have 
the form u{t) = exp(/j,t)0(t) where (j){t) is a periodic func- 
tion with period ^-klo^Jv |32]- One can readily see that 
the solutions u{t) are unstable if Re(/i) > 0. We per- 
formed a perturbative stability analysis to determine the 
stability boundary for the solutions [32]. We find that 
the instability appears when the condition 



2^0 1 



(19) 



is fulfilled. This leading order expression is precisely 
twice the excitation frequency of Eq. (pi). Finally, we 
check our result by solving the system of coupled non- 
linear equations (5a) and (5b) in real time. In Fig. [51 



we plot for each modulation frequency v the maximum 
number of photons after the oscillation has stabilized, 
starting with very small initial (a, /3). Thus when the 
frequency of modulation is chosen correctly, there will be 



a substantial photon flux at the cavity face, even below 
threshold. We should also note that there is a broadening 
in resonance modulation frequencies as the critical point 
is approached (see Fig. Isl), and this is consistent with 
the behavior of the imaginary part of the eigenmode in 

Fig.m 





(a) 



(b) 




FIG. 5: (a) The maximum number of photons JQp and (b) 
the maximum value for Re(/3), as functions of the coupling 
strength A and the modulation frequency v. The peaks occur 
at i^ = 2Re(tJe2,) where uje^. is given by Eq. (|9|. In both figures 
the red curves show the instability condition Eq. ( 19 1, (c) The 
oscillations of Iie{P)/N in time are shown for A = O.SAc where 
the modulation frequency is chosen to be equal to resonance 
frequency given by z^ = 1.2a;o. For all figures here, e = 1/50 
and system parameters are same as in Fig. [2] 



that the coupled BEC-cavity system can be mapped into 
an open-system realization of the Dicke model with a 
symmetry-breaking field. 

The zero-field phase transition is driven by softening of 
polaritonic excitations, which provides access to the in- 
ternal dynamics of the coupled system close to the crit- 
icality. We discuss a number of photodetection-based 
techniques for the photons that leak out of the cavity 
mirrors and relate it to intra-cavity critical dynamics. 

We find that the intra-cavity photon number, which 
can be measured by a photodetector outside the cavity, 
displays a different scaling law for the open system than 
the closed system with perfectly reflecting mirrors. In the 
latter case, the dynamics conserves total excitation num- 
ber and the photon number is calculated in the ground 
state. 

We next discussed the second order correlation func- 

(2) 

tion for photons gls (t) in the steady state, via coinci- 
dence measurement of photons leaking out of the cavity 
walls. We show that gss {t) displays damped oscillations 
with a period that increases progressively as the critical 
point is approached, signalling the critical slowing down 
of the coupled intra-cavity dynamics. We relate the spec- 
tral content of the time-series of gis (t) to the complex 
frequency of the softening polaritonic mode. We also 
show that trap misalignment can lead to a background 
coherent cavity field and that its signature is a charac- 

(2) 

teristic beating pattern in gss {t). 

Finally, we discussed a modulation scheme that di- 
rectly captures the softening of the polaritonic mode. 
This is done by introducing a parametric instability in 
the system through the periodic modulation of the drive 
Rabi frequency in time. We show that modulation at 
twice the polaritonic mode frequency results in a measur- 
able photon flux outside the cavity below the threshold. 



V. CONCLUSION 

In this work, we discussed the excitations of an opti- 
cally driven atomic condensate coupled to a single mode 
of a high-finesse cavity that displays a self-organization 
transition as a function of the driving strength. Tak- 
ing into account the finite-size of the system, we showed 
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